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Summary. The propagation of solitary waves in a Bragg grating formed by an 
array of thin nanostructured dielectric films is considered. A system of equations of 
Maxwell-Duffing type and describing forward- and backward-propagating waves in 
such a grating, is derived. Exact solitary wave solutions are found, analyzed, and 
compared with the results of direct numerical simulations. 



1.1 Introduction 

The last decade has been a period of rapid progress in the field of photonic 
crystals [1]-[10] . In particular, the one-dimensional case of a resonant Bragg 
grating [l]-[5] or a resonantly absorbing Bragg reflector (RABR) [6]-[8] has 
been studied extensively. In the simplest case a resonant Bragg grating consists 
of a linear homogeneous dielectric medium containing an array of thin films 
with resonant atoms or molecules. The thickness of each film is much less 
than the wavelength of the electromagnetic wave propagating through such a 
structure. The interaction of ultra-short pulses and films embedded with two- 
level atoms has been studied by Mantsyzov et al. [l]-[5] in the framework of the 
two- wave reduced Maxwell-Bloch model and by Kozhekin et al. [6]- [8]. This 
work demonstrated the existence of the 2-7r-pulse of self-induced transparency 
in such structures [1, 4, 6]. It was also found [8] that bright as well as dark 
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solitons can exist in the prohibited spectral gap, and that bright solitons can 
have arbitrary pulse area. 

If the density of two-level atoms is very high, then the near-dipole-dipole 
interaction is noticeable and should be accounted for in the mathematical 
model. The effect of dipole-dipole interaction on the existence of gap solitons in 
a resonant Bragg grating was studied in [10] ; details can be found in [9] . Recent 
numerical simulations have yielded unusual solutions known as zoomerons. [19] 
The optical zoomeron was discovered and investigated recently [20] in the 
context of the resonant Bragg grating. A zoomeron is a localized pulse similar 
to an optical soliton, except that its velocity oscillates about some mean value. 
These works also contain careful construction of the underlying mathematical 
model, which is derived from first principles. 

Recent advances in nanofabrication have allowed the creation of nanocom- 
posite materials, which have the ability to sustain nonlinear plasmonic oscilla- 
tions. These materials have metallic nanoparticlcs embedded in them [21]-[23]. 
In this paper we consider a dielectric material into which thin films containing 
metallic nanoparticles have been inserted. These thin films are spaced period- 
ically along the length of the dielectric so that the Bragg prohibited spectral 
gap is centered at the plasmonic resonance frequency of the nanoparticles. 
We derive governing equations for slowly-varying envelopes of two counter- 
propagating electromagnetic waves and of the plasmonic oscillation-induced 
medium polarization. We find that this system of equations has the form of the 
two-wave Maxwell-Duffing model. We find exact solutions of this system and 
demonstrate that, in contrast to conventional 27r-pulses, they have nonlinear 
phase. We show that the stability of these solutions is sensitive to perturba- 
tion of this phase. We also study the collisions of these pulses and find that 
the outcomes of such collisions are highly dependent on relative phase. 

1.2 Basic Equations 

We consider a grating formed by an array of thin films which are embedded 
in a linear dielectric medium. In our derivation of the governing equations we 
follow [l]-[8], wherein Bragg resonance arises if the distance between successive 
films is a = (A/2) m, m = 1, 2, 3, .... To obtain governing equations we apply 
the transfer-operator approach, presented below. 

1.2.1 Transfer-operator approach 

Let us consider the ultra-short optical pulse propagation along the X-direction 
of the periodic array of thin films, which are placed at points . . . , x n -i, x ni x n+ \ 
(Fig.l). The medium between films has dielectric permittivity e. Hereafter, for 
the sake of definitcness, we consider a TE-wave whose electric field component 
is parallel to the layers. All results can be generalized easily for the case of 
TM-polarized waves. 
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It is suitable to represent the electric and magnetic strengths E, if, and 
the polarization of the two-level atoms ensemble P in the form of Fourier 
integrals 

oo 

E(x,z,t) = (27r)~ 2 J cxp[-iLut + if3z}E(x, [3,Lu)dt dz, 

— oo 
oo 

H(x, z, t) = (27r)~ 2 J exp[-iuut + i(3z]H(x, (i, u)dt dz, 

— oo 
oo 

P(x n ,z,t) = (27r)~ 2 J cxp[—iLut + if3z}P(x n ,f3,Lu)dt dz. 

— oo 

Outside the films the Fourier components of the vectors E(x,f3,u>) and 
H{x,j3,Lu) are defined by the Maxwell equations. At points x n these values 
are defined from continuity conditions. Thus, the TE-wave propagation can 
be described by the following system. 

g + (^-/3 2 )£ = 0, (1.1) 

H x = -((3/k)E, H z = -(i/k)dE/dx, E V = E, 

with boundary conditions [11, 12] 

E(x n -0) = E(x n + 0), H z (x n + 0)- H z (x n -0) = 4iTrkP y (x n ,p,w), (1.2) 

where k = lu/c. The solutions of equation (1.1) in the intervals x n < x < x n+ \ 
can be written as 

E(x, /3, lu) = A n (/3, lu) exp[iq(x - x n )\ + B n (j3, lu) exp[-iq(x - x n )], 
H z (x, P, lu) = qk^ 1 {A n (/3, lu) exp[iq(x - x n )] - B n {(3, lu) exp[-iq(x - x n )}} , 



where q = \J k 2 e — [3 2 . Hence, the amplitudes A n and B n completely deter- 
mine the electromagnetic field in a RABR. Let us consider the point x n . The 
electric field at x — x n — 5 (5 < < a) is defined by amplitudes and B^ , 
and the field at x — x n + S is defined by A^ and B^ . Continuity conditions 
(1.1) result in the following relations among these amplitudes 

AW +^ =A f) 
AW -BW=AM -sW +4 7rffcV 1 P5, n , 

where Ps, n = Ps{A^ + B^) is the surface polarization of a thin film at 
point x n , which is induced by the electrical field inside the film. Thus we find 
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A { n R) =A^+27rik 2 q- 1 P s , n , B™ = B n L ^ — 2mk 2 q~ 1 Ps, n - (1-3) 

Taking into account the strength of the electric field outside the films, we 
write 

A n L +i = 4*>exp(iga), B n % = B n R ^ ex V (-iqa). (1.4) 

If the vectors tp n L>> = (A n L \ B n L ^) and ip n R '^ — (An®, B n R ^) are introduced, 
then the relations (1.3) can be represented as 

where U n is the transfer operator of vector ip n L ^ through the film located at 
point x n . In the general case U n is a nonlinear operator. The relations (1.4) 
are represented in the vectorial form 

where the linear operator V n transfers the vector ip n R ^ between adjacent thin 
films and is represented by the diagonal matrix 

f> _ ( exp(iqa) \ 
Vn ~ \ exp(-iga) J ■ 

In this manner we define the nonlinear transfer-operator of the vector 
through an elementary cell of RABR: 

^n%=VnU n ^=f n ^n L) - (1-5) 

The transfer-operator approach is frequently used in models of one-dimensional 
photonic crystals of linear media, e.g. distributed feedback structures [13]. 

In (1.5) the upper index can be omitted, and the equation can be rewritten 
as the following recurrence relations 

A n+ i = A„, exp(iga) + 2mk 2 q~ 1 P s , n exp(iqa), (1.6) 
B n+ i = B n cxp(-iqa) - 2Trik 2 q^ 1 P s ^ n cxp(-iqa) (1.7) 

These recurrence relations are exact, as no approximations (e.g. slowly-varying 
envelope of electromagnetic pulses approximation, the long-wave approxima- 
tion) have been employed. Furthermore, the surface polarization of a thin film 
could be calculated via different suitable models. Here we follow the works by 
B. Mantsyzov at al. [l]-[5], and A. Kozhekin, G. Kurizki, et al. [6]- [9], where 
the two-level atom model has been used. 

1.2.2 Linear response approximation 

To demonstrate that the RABR is a true gap medium, it is suitable to obtain 
the electromagnetic wave spectrum through a linear response approximation. 
In the general case we can use the following expression for polarization: 
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P s , n =xH(4 fl) + B( n H) )- (1-8) 
Substitution of this formula in (1.6), (1.7) yields 

A n+ i = (1 + ip)A n exp(iqa) + ipB n cxp(iqa), (1.9) 
B n+ i = (1 - ip)B n exp(-iqa) - ipA n exp(-iga). (1-10) 

Here p = p(w) = 2-Kk 2 q~ 1 x{ijj) = 27ra>c~ 1 £~ 1/ ' 2 x(u>). We employ an ansatz 
in which the wave is a collective motion of the electrical field in the grating. 
Hence 

A n — Acxp(ikna), B n = B cxp(ikna). (1-H) 

Refs. (1.9), (1.10) show that the wave amplitudes A and B satisfy the following 
linear system of equations: 

Aexp(iKa) = (1 + ip)Acxp(iqa) + ipB exp(iqa), (1-12) 
B exp(iKa) = (1 — ip)B cxp(— iqa) — ipAcxp(-iqa). (1-13) 

A nontrivial solution of this system exists if and only if the determinant is 
equal to zero, i.e. 

^^f (1 + ip) cxp(iqa) — cxp(iKa) ipexp(iqa) 

y —ipcxY>{—iqa) (1 — ip) cxp(— iqa) — exp(iKa) 

(1.14) 

If we define 

Z = cxj>(iKa), 

G = (1 + ip) exp(iqa) — (cos qa — psinqa) + i(pcosqa + sinqa), 

then equation (1.14) can be rewritten as the following equation in Z: 

Z 2 - (G + G*)Z+l = 0. 

This equation has solutions 

Z± = Re G ± iy/l - (Re G) 2 . 

If Re G < 1, then Re G ± iy/l - (Re G) 2 = cos Ka + i sin Ka. Hence the 
wave numbers K± are real-valued and satisfy the transcendental equation 

cos Ka — cos qa — psinqa. (1-15) 

If Re G > 1, then the roots of equation (1.14) are real. In this case, the 
wave numbers K± are pure imaginary. The condition Re G > 1 defines the 
frequencies of the forbidden zone. The waves with these frequencies cannot 
propagate in the grating. The boundaries of this forbidden zone are defined 
by 
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Fig. 1.1. An example of band structure, if inclusions are linear oscillators. 



cos qa — p sin qa = 1 . (1-16) 

The model of the resonant system containing the thin films defines the 
explicit form of the function p = p(cj). The form of p{uo) determines the 
dispersion relation (1.15). It should be noted that this dispersion relation 
ensures a series of gaps in the electromagnetic wave spectrum. Figure (1.1) 
represents an example of such a band structure when the inclusions are linear 
oscillators. 

1.2.3 Long-wave and weak nonlinearity approximations 

By using the long-wave approximation we can transform the exact equations 
(1.6), (1.7) into differential equations. To do so we introduce the field variables 

A( x ) — X! A n 5(x-x n ), B{x) = y~]B n 6(x-x n ), P(x) = ^ Ps, n 5{x-x n ), 

n n n 

Using the integral representation for a Dirac delta-function, 

oo 

S(x) = (2n)^ 1 J exp(ikx)dk, 

— oo 



we produce the following expression 
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oo 

A(x) = (27r) _1 ^A„ / exp[ik(x - x n )]dk 

OC 

= (2.)- ] eM*kx)Y,A n eM-*kx n)dk . 

-oo 

It follows that the Fourier transform of A(x) is then 

A(k) — A n cxp(—ikx n ) — A n exp(-ikan) , 

n n 

and the Fourier transforms of B{x) and P(x) are 

B(k) = B n exp(— ikan), P{k) = Ps,n exp(— ifccm), 

n n 

respectively. 

The form of these spatial Fourier components ensures periodicity in k. For 
example, 

A(k) = ^ A n cxp(-ikan) exp(±27rm) 

n 

= ^ A n exp(—ikan ± 2-Kin) (1-17) 

n 

= A n exp[-mn(fc ± 27r/a)] = A(k ± 27r/a). 

n 

^From recurrence equations (1.6), (1.7) we have 

A(k) exp(ika) — A(k) exj>(iqa) + inP(k) exp(iga), 
£>(fc) exp(i/ca) = B(k) cxp(-iqa) — inP(k) exp(— iqa), 

or 

[exp{ia(fc - q)} - 1] A(k) = inP{k), (1.18) 

[cxp{m(fc + q)} - 1] B{k) = -inP{k). (1.19) 

The linear approximation P(k) — x(cj)[j4(fc) + B(k)] for polarization 
again leads us to the linear dispersion law (1.14). The system of equations 
(1.18), (1.19) is equivalent to the discrete equations (1.6), (1.7). At this point 
we have made no approximations other than the assumption on the thin films' 
width. 

If the thin film array were absent, then the dispersion relation would be 

cosfco = cos qa. 
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In this case the wave with amplitude A(k) has k = q, indicating propagation 
to the right, while the wave with amplitude B{k) has k = —q, indicating 
propagation in the opposite direction. Let us suppose that the polarization of 
the thin film array produces little change in wave vectors, i.e., for the right- 
propagating wave the wave vector is k = q + 5k, and for the opposite wave 
the wave vector is k = — q + 8k. If we set the value of q near one of the Bragg 
resonances, say q = 2n/a + 5q, where Sq <C 2ir/a, then equations (1.6), (1.7) 
take the form 

[exp(ia5k) - 1] A(q + Sk) = inP(q + 5k), 
[exp(iaSk) - 1] B(-q + Sk) = -inP(-q + 5k). 

By virtue of the periodicity conditions (1.17) these equations can be re- written 

as 

[exp(ia5k) - 1] A(5q + Sk) = iK,P(Sq + Sk), 
[cxp(iaSk) - 1] B(-5q + Sk) = -inP(-Sq + Sk). 

After the change of variables Sk = Sk ± Sq we have 

exp{ia(Sk - Sq)} - lj A(Sk) = inP(Sk), (1.20) 

exp{m(<5fc + Sq)} - lj B(5k) = -inP{S~k). (1.21) 

The long- wave approximation means that nonzero values of the spatial Fourier 
amplitudes are located near zero value of argument. Suppose that aSk is small 
enough that e M$k+Sq) ~ ^ + ia(Sk + Sq). Then we replace equations (1.20), 
(1.21) with the approximate equations 

ia(Sk - 5q)A{5k) = i K P(Sk), (1.22) 

ia(Sk + Sq)B(Sk) = -inP(Sk). (1.23) 

If now we return to the spatial variable, equations (1.22), (1.23) lead us to 
the equations of coupled- wave theory 

dA 

— = i8qA(x) + ina^Pix), (1.24) 
3B 

— = -iSqB(x) - ina^Pix). (1.25) 

In these equations the fields A(x), B(x), P(x) and the parameters Sq, n 
are functions of frequency to. To obtain the final system of equations in the 
spatial and time variables, we perform an inverse Fourier transformation. We 
now assume that the envelopes of the electromagnetic waves vary slowly in 
time [14]. This approximation simplifies the system of coupled wave equations 
under consideration. 
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1.2.4 Slowly-varying envelope approximation 

In the slowly-varying envelope approximation, we assume that the approxi- 
mated fields are inverse Fourier transforms of narrow wave packets, i.e. they 
are quasiharmonic waves [14]. For example, the quasiharmonic wave form of 
the electric field is 

E(x,t) = £ (x,t) exp[— iojot], 

where ujq is carrier wave frequency. The electric field E(x, t) and the Fourier 
components of the envelope £(x, t) of the pulse are related by 

oo 

E(x,t) — (27r) _1 J E(x,u>) exp(—iojt)dbj 

— oo 
oo 

= (27r) _1 J £(x, uj) exp[— i(oj + ujo)t}duj 

— oo 
oo 

= (27r) _1 J £(x, uj — w ) exp(— iujt)duj, 

— oo 

where the function E(x, ui) is nonzero if u> e (oj n — Au), uj + Au), with Aw <C 
loq- Hence, E(x,uj + ujq) — £(x,u>). Thus, if we have some relation for E(x,uj), 
then the analogous relation for £(x, uj) can be found by shifting uj ^> uj + uj 
in all functions of uj. 
Let 

A(x,t) — A{x,t) exp(— iuJot), B(x,t) — B(x, t) cxp(— iuj^t), P(x,t) — V(x, t) exp(— iu>ot). 

^From Eqs. (1-24), (1.25) it follows that the Fourier components A(x, uj),B(x, uj), 
and V(x,uj) satisfy 

dA 

-g—(x,uj)=i5q(ujQ+uj)A(x,uj)+iK(uj +uj)a V(x,uj), (1-26) 
— (x, uj) = —iSq(u) + uj)B(x, uj) — zk(u + uj)a^ lr P(x, uj) (1-27) 
Since A, B and V are nonzero for uj <C uj , one can use the expansions: 

Sq(ui + uj) ~ <7o — 2n/a + qiuo + q 2 uj 2 /2, k(w + w)a _1 m K n , (1.28) 

where q n = d n q/du> n at u> = uj - In particular, q^ 1 = v g is the group velocity, 
and <72 takes into account the group-velocity dispersion. 

Considering expansions (1.28), we have the following description of the 
evolution of slowly-varying envelopes. 
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j (l;-^l) E+ lS E -^ 8 = +w <L30 > 

where Aq = qo — 2ir / a. The next step requires a choice of the model for 
the thin films' medium. Possibilities include anharmonic oscillators, two- or 
three-level atoms, excitons of molecular chains, nano-particles, quantum dots, 
and more. First we consider the two-level atom model. 



1.2.5 Example: thin films containing two-level atoms 

Here we employ the approach developed above to derive an already known 
system of equations [1]. The model assumes that each thin film contains two- 
level atoms. The state of a two-level atom is described by a density matrix p. 
The matrix element pi2 describes the transition between the ground state |2) 
and excited state |1). P22 and p\\ represent the populations of theses states. 
Evolution of the two-level atom is governed by the Bloch equations [15]. 

d 

l hgj_Pi2 = hAjjpu - d 12 (p 2 2 - Pu)A in , (1.31) 

d 

tti—^p-n - pn) = 2 {d 12 p2iA in - d 2 ipi 2 A* n ) . (1.32) 

In these equations A in is the electric field interacting with a two-level atom. 
In the problem under consideration, A in = A + B, and 

T . „ 2iruj( 3 n at d 1 2 , . 

K V = -. — 7 — P12 • 

cn(w ) 

Here the cornerstone brackets denote summation over all atoms within a fre- 
quency detuning of Auj from the center of the inhomogeneity broadening line, 
n(u>o) is the refractive index of the medium containing the array of thin films, 
and n at is the effective density of the resonant atoms in the films. n at is de- 
fined by n at — N at (£f/a), where N at is the bulk density of atoms, if is the 
film width, and a is the lattice spacing. 

We suppose the group-velocity dispersion is of no importance. The result- 
ing equations are the two-wave reduced Maxwell-Bloch equations. We intro- 
duce the normalized variables 

ei = todi 2 A/h, e 2 = t a di 2 B/h, x = (v g t a , T = t/to- 

The normalized two-wave reduced Maxwell-Bloch equations take the following 
form: 

( d d\ 
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( d d\ 

i ( q£ ~ -7^ J e2 - Se 2 = +7<Pi2) , (1-34) 
d 

% 'q~P^ = Api2 - ne in , (1.35) 
d 

—n = -4 Im (pi2e in ), (1.36) 

where 7 = t v g /L a , 5 = t v g Aq , L a = (cn(u> Q )h) / (2Tru> t a n a t\d 12 \' 2 ) is the 
resonant absorption length, and A = Auto is the normalized frequency de- 
tuning. 

We define n = P22 — P11, &in = e i + ^2 and introduce yet another change 
of variables: 

e»„ = ei + e 2 = / s exp(i<5r), ei - e 2 = / a cxp(«<$T), p 12 = rcxp(i5r). 
The system of equations (1.33) - (1.36) can be rewritten as 

f + f =0, (1.37) 

^ + |^ = 2z 7 (r), (1.38) 

i^-r = (^ + J)r-n/ a , (1.39) 

|-n = -4 Im (r/ 8 *). (1.40) 

/From (1.37) it follows that 

dfa = _dh 

dC, dr ' 
which allows to rewrite (1.37)-(1.40) in the form 
d 2 f s d 2 .f s . /Or 



oe ^=- 2i n^/' (L41) 

i— r = (4 + S)r - nf s , (1.42) 

— n = -4 Im (r/ a *). (1.43) 

If we assume that the inhomogeneous broadening is absent, i.e. if the 
hypothesis of a sharp atomic resonant transition is true, then 6 + A = 0, and 
this system reduces to the Sine-Gordon equation [1]. Reference [4] presents the 
steady-state solution of (1.41)-(1.43) with inhomogeneous broadening taken 
into account. 
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1.2.6 Thin films containing metallic nanoparticles 

It was shown above that counter-propagating electric field waves A and B 
in the slowly-varying envelope approximation satisfy the following system of 
equations: 

where Z\<?o = <Zo — 2n/a is the mismatch between the carrier wavenumber and 
the Bragg resonant wavenumber. To describe the evolution of material polar- 
ization in the slowly-varying amplitude approximation, we must model the 
thin films' response to an external light field. Previous work has considered 
various mechanisms as sources of the dielectric properties of mctamaterials. 
In the simplest case, dielectric properties can be attributed to plasmonic oscil- 
lations, which arc modeled by Lorcntz oscillators. Magnetic properties can be 
described by the equations of a system of LC-circuits [27, 28, 29, 30, 31]. The 
simplest generalizations of this model include anharmonicity of plasmonic 
oscillations [16, 32] or the addition of a nonlinear capacitor into each LC- 
circuit [33]. In this paper we consider an array of non-magnetic thin films 
containing metallic nanoparticles, which have cubic nonlinear response to ex- 
ternal fields [21, 22, 16]. 

The macroscopic polarization P is governed by the equation 

d 2 P , dP , ojI 

-W + ^ P + r «lH +KP = 4^ 

where uj p is plasma frequency and iOd is dimension quantization frequency 
for nanoparticles. Losses of the plasmonic oscillations are taken into account 
by the parameter P a . It is assumed that the duration of the electromagnetic 
pulse is small enough that dissipation effects can be neglected. If the anhar- 
monic parameter n is equal to zero, then we have the famous Lorentz model 
for describing electromagnetic wave propagation and refraction in metamate- 
rials [35, 36, 27, 28, 29, 30, 31]. 

Starting from the slowly-varying envelope approximation, standard ma- 
nipulation leads to 

+ ("d - "o)V + ^\P\ 2 V - --p^£ int (x,t). (1.46) 
at Iujq airujQ 

Terms varying rapidly in time, which are proportional to exp(±3zwo^); are 
neglected. In this equation £i nt is the electric field interacting with metallic 
nanoparticles. In the problem under consideration we have £i nt = A + B . 

Due to the limitations of nanofabrication, the sizes and shapes of nanopar- 
ticles are not uniform. In practice, deviation from a perfectly spherical shape 
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has a much larger impact on a nanoparticle's resonance frequency than does 
variation in diameter. This causes a broadening of the resonance line. The 
broadened spectrum is characterized by a probability density function g(Aui) 
of deviations Aw from some mean value to res . When computing the total po- 
larization, all resonance frequencies must be taken into account. 

The contributions of the various resonance frequencies are weighted ac- 
cording to the probability density function g(Acj); the weighted average is 
denoted by (V) in equations (1.44), (1.45). In what follows, n{u>o) denotes 
the refractive index of the medium containing the array of thin films, and 
n np is the effective density of the resonant nanoparticles in films. As in the 
model of films containing two-level atoms, the effective density is equal to 
n np — N np (£f/a), where N np is the bulk density of nanoparticles, if is the 
width of a film, and a is the lattice spacing. 

We study a medium-light interaction in which resonance is the dominant 
phenomenon. As such, the length of the sample is smaller than the charac- 
teristic dispersion length. In this case the temporal second derivative terms 
in equations (1.44,1.45) can be omitted. The resulting equations are the two- 
wave Maxwell-Duffing equations. They can be rewritten in dimensionless form 
using the following rescaling: 

ei = A/Ao, 
e 2 = B/A , 

p = (4iruo/[^/euj p Ao})'P, 
C = {u p /2c)x, 

T = t/to- 

Here to = 2^/s/lo p , while Aq is a characteristic amplitude of counter- 
propagating fields. In dimensionless form, the two- wave Maxwell-Duffing equa- 
tions read 



d_ d_ 
d__d_ 



ei + Sei = -(p), 

e 2 -Se 2 = +<p), (1.47) 



.dp l2 i \ 

+ Ap + fi\p\ p = -(ei + e 2 ), 

where /i = (3Kv / e/woWp)(v / £Wp/47rwo) 2 ^4o ^ s a 
anharmonicity, 5 = 2Aqo(c/uj p ) is the dimensionless mismatch coefficient, 
A = 2y/e(uJd — uoq)/lo p is the dimensionless detuning of a nanoparticle's reso- 
nance frequency from the field's carrier frequency. 

In a coordinate system rotating with angular frequency S, 

ei=/ie* 5r , e 2 = he UST ,p = qe itT , 



equations (1.47) become 
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■ fd__9_\ 
' ' d( dr) 



h = +<«>, (1-48) 



i^ + (4-J)g + p|g| 2 g = -(/1 + /2). 

Further simplification of the system (1.48) can be achieved by introducing 
new variables 

fs = ~(fl + f2), fa = fl — /2i 

which allow decoupling of one equation from the system of three equations. In 
these new variables the polarization q is coupled with only one field variable. 
Simple transformations give 



(q), (1.49) 
(q), (1.50) 



d 2 f a 


d 2 fa 


„ 9 


d( 2 


dr 2 


d( 


d 2 f s 
d( 2 


d 2 .f s 
8t 2 


n 9 



. dq 



(A-S)q + v\q\ 2 q = fs- (1.51) 
As one can see, we have a coupled system of equations for f s and q. 



l dr 



1.3 Solitary Wave Solutions 

We consider localized solitary wave solutions of (1.49) in the limit of narrow 
spectral line Aujg/Auj s <C 1, where Alo s and Auj g are spectral widths of a signal 
and spectral line g(Au>). In this case the spectral line can be represented as 
Dirac ^-function: g(Aui) = 5(Aui). Equations (1.49), (1.51) can then be re- 
written as follows: 

l^- + (A-5)q + n\q\ 2 q = f, (1.53) 

Scaling analysis of this system shows that solitary wave solutions can be rep- 
resented as 



1 / 2v 2 



3/4 



f s = f F n{ r ] ) = -^) Fair,), 
1 / 2v 2 \ 1/4 



V = (C - vt) 



1 — V 



2 ' 
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Here v is velocity of the solitary wave, r\ is a scale-invariant parameter in 
a coordinate system moving with the solitary wave, and functions Fq, Qq 
satisfy the following system of equations: 

F" = -iQ', (1.55) 
-iQ' + HQ + \Q\ 2 Q = F. (1.56) 

The only dimensionless parameter which remains in the system is 

n=(A-5) ] f^f, (1.57) 

which characterizes the deviation of carrier frequency from the plasmonic 
frequency u> p and the Bragg resonance frequency to Br- 

The first equation implies that F' = —iQ + constant. We seek a solitary- 
wave solution, so we assume that F, Q, and their derivatives decay to zero as 
\r)\ — > oo. Hence the constant is zero, and 

iF' = Q 

-iQ' + t2Q + \Q\ 2 Q = F. (1.58) 
The system of ordinary differential equations (1.58) has integral of motion 
\Q\ 2 - \F\ 2 = constant. 
For solutions decaying as \rj\ — > oo, the constant is equal to zero, and 

|Q| 2 HF| 2 . 

This allows the following parametrization of solutions: 

F(rj) = R( V )e^) ; Q^) = R( v y^) ; 

where R, </>, and ip are real-valued functions satisfying 

R' = -Rsm(<t>-i>), 

<f>' = -cos(0- V), (1-59) 
ip' + Q + R 2 = cos((f>- 

If we set $ = 4> — tjj, then we have 

$' - Q- R 2 = -2cos<£, 

R' = -Rsm$. (1.60) 

Taking into account second equation of (1.60), the first equation can be re- 
written as 

d 

R— (cos<P) - Q - R 2 = -2cos#. 
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If we set y = cos <P, then we have 

R% + 2y = R 2 + n. (1.61) 
art 

The solutions of equation (1.61) have the form 

y = \R 2 + j + cR- 2 , 

where c is free. Since y = cos<P, the right-hand side must remain between 
-1 and 1. As we expect R — > as \rj\ — > oo, c must be zero. We have the 
conservation law 

cos<P = -R 2 + — . (1.62) 
4 2 

Substituting (1.62) into the second equation of (1.60) and subsequent integra- 
tion gives the following expression for R: 

R 2_ 2(4 -n 2 ) 



Q + 2 cosh { VI - 7F(77 - no) + \ In (j^pr) } ' 

The right-hand side is positive real- valued for all -q if and only if — 2 < fi < 2. 
\/4 — i? 2 appears so often in what follows that we set @ = \J\— fi 2 . R 2 then 
has the form 

R2 = ^ + 2cosh{/3(r/-?7')}' (L63) 
where we have combined arbitrary constant r] and the logarithm into a single 
constant (77') in the argument of the hyperbolic cosine. Using the conservation 
law (1.62) we obtain an expression for <£: 

<2> = 2 arctan [ ~ tanh j ^(77 - 77') 1 ) . (1.64) 



Now we integrate (j)' = — cos and find 

^/ a /2- , r 1 „, ,, 

<p=— — (77 — ?7 ) — arctan — - — tanh < -p(?7 — 77 ) 



2 w " V /3 12 

Finally, we determine ip: 

ip = 4> — <p 

= —-^{V 'H ) ^ arctan I — — — tanh < — 7/1 ' 

This pulse exists only if value of the parameter J? are inside the interval 
— 2 < fl < 2. The maximal value of the amplitude of this solitary solution is 

A = y/2(2-n). 
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The phases 4> and i 1 are nonlinear. Their behavior is asymptotically linear as 
77 — > ±00. If SI = 0, then the limiting values of the phases satisfy 

|0(oo) - <f>(-oo)\ = it/2, 
|V>(oo) - V(-oo)| = 3tt/2. 



1.4 Energy Partition 

The total energy of the solitary wave is distributed among co/contr-propagating 
fields and medium polarization. Here we study energy partition between all 
these components. Using equations (1.49), (1.50) and conditions as — ► 00, 
one can show that 

fain) = -\f s {n) = ~F{ n ). (1.65) 

We are interested in the energies of the dimensionless fields f\ , f 2 , and polar- 
ization q. 

h = \(fa-fs) = -^ (l + 1 ) F 

/„ = -!(/. + /„) = § (1-i)f 
q = qoQ 

To find energies of forward- and backward-propagating waves we need to cal- 
culate an energy of a solitary wave: 



+00 

E R = J \F\ 2 ( V )dr, = J " ' ~ 



+00 



R 2 (r])df] = 8 arctan \j * - ^ . (1.66) 



Finally we have the energies 



4 \v 

E /2 = |(i-l) 2 ^, (1.68) 

E q = q 2 E R . (1.69) 

Ratios of energies in different fields as well as polarization have the following 
form 
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Eh 




E h 


i(i+«) 


E q 


2(l-vY 


Ef 2 


l(l-v) 


E q 


2(l+«)' 



(1.70) 
(1.71) 
(1.72) 

Therefore, energy partitioning is determined by only one parameter v which 
is dimensionless combination of main system parameters. 



1.5 Numerical simulation 

The shape and phase of the incident pulse are controllable in real experimental 
situation. To model pulse dynamics in the Bragg grating it is natural to con- 
sider asymptotic mixed initial-boundary value problem for equations (1.48). 
We define initial conditione as 

«(C,t)->0, /i(C,t)->0, /2(C.r)-0, r^-oo (1.73) 

with no incident field at the left edge of the sample and with incident field at 
the left edge defined as follows: 



fi (— 10, t) = wexp(it 



w = 3.5 exp 



1 /t-3.0 n 



2 V 1-5 
6 = arctan (tanh [1.5 (t - 3.0)]) . 



(1.74) 



In our case the spatial simulation domain was chosen as [—10,40]. Parameters 
A — S and /j, were chosen 

A- 6 = 0, n=l. (1.75) 

As one can see we gave the initial pulse the same configuration of phase in 
topological sense as in solitary wave solution. This point is important, because 
otherwise phase difference cannot relax to the symmetry of the stationary wave 
which is revealed in (1.64). As a result without right "topological charge" 
solution will be unstable. 

For the second field boundary condition was as follows 

/ 2 (40,r) = 0. (1.76) 

In the first experiment we injected a pulse relatively close to solitary wave 
solution. The results are shown in Figs. 1.2-1.4. The amplitude phase differ- 
ence slightly differed as well as pulse shape which was Gaussian. On the first 
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T 

Fig. 1.2. Propagation of pulse. The first experiment. Mapping of the |ei(£,T)| 
surface. 




T 

Fig. 1.3. Propagation of pulse. The first experiment. Mapping of the |e2(C, t)| 
surface. 



stage (t < 7) of evolution we observed fast excess energy damping in radia- 
tion of quasi-linear waves in both directions and relaxation to solution roughly 
close to stationary one. Then we had some stage of pulse shape refinement 
(7 < t < 30) with consequent propagation of the solution very close to (1.63). 
One can compare refined pulse shape with stationary solution in Fig. 1.5. 
During second experiment we used pulse of the lower amplitude 



fi (— 10, r) = wexp(i6 



w = 2.0 exp 



1 /t-3.0 

2 V 1-5 



6 = arctan (tanh [1.5 (r - 3.0)]) 



(1.77) 
(1.78) 




t 

Fig. 1.4. Propagation of pulse. The first experiment. Mapping of the |g(£,r)| sur- 
face. 
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Fig. 1.5. The first experiment. Absolute value of polarization at different moments 
of time. Stationary pulse shape almost coincides with exact solution. 



Results are represented in Figs. 1.6-1.8. Clearly observable oscillations of 
pulse magnitude can be associated with modes localized on pulse. This very 
instructive result could be a fruitful ground for further investigation and ex- 
periments. 
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Fig. 1.6. Propagation of pulse. The second experiment. Mapping of the \ei((,r) 
surface. 
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Fig. 1.7. Propagation of pulse. The second experiment. Mapping of the |e2 (Ci 7 ~) 
surface. 



1.6 Conclusion 

We derived equations describing optical pulse evolution in Bragg gratings 
with thin films containing active dopants. In particular we investigated case 
of thin films containing metallic nanoparticles. We showed that corresponding 
system have one parameter. This parameter contains information on Bragg 
and plasmonic frequencies and general pulse characteristics. These equations 
have solitary wave solutions describing bound state of two contraindicating 
waves with medium polarization. 
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Fig. 1.8. Propagation of pulse. The second experiment. Mapping of the |p(Ci r )l 
surface. 
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